Partial synchronisation of stochastic oscillators through hydro dynamic coupling 
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Holographic optical tweezers are used to construct a static bistable optical potential energy land- 
scape where a Brownian particle experiences restoring forces from two nearby optical traps and 
undergoes thermally activated transitions between the two energy minima. Hydrodynamic cou- 
pling between two such systems results in their partial synchronisation. This is interpreted as an 
emergence of higher mobility pathways, along which it is easier to overcome barriers to structural 
rearrangement. 
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Activated dynamics is ubiquitous in physics, chemistry 
and biology pi, f2|. Escape from a stable state occurs via 
transition events by which a system crosses an energy 
barrier, falling into a neighboring state. On the small 
length scales of colloidal particles and biological macro- 
molecules, it is the surrounding fluid which provides the 
required fluctuations to overcome a barrier and, at the 
same time, the viscosity needed to relax the system in the 
second minimum by dissipating energy. With the addi- 
tion of a spatial potential landscape consisting of two 
or more minima, such an environment allows a colloidal 
particle to undergo hops between the minima, provided 
that any energetic maximum is surmountable by ther- 
mal fluctuations [3]-[5]. Starting from a Fokker-Planck 
formulation of Brownian motion in bistable potentials, 
Kramers [1 was able to obtain an analytical expression 
for the hopping rates that depends solely on the local fea- 
tures of the underlying potential. McCann et al. [3 used 
two separate HeNe lasers to construct a bistable optical 
landscape where a 600 nm Silica sphere was confined in a 
solvent fluid. The sphere hopped back and forth between 
the traps, providing the ideal representation of Kramers' 
ideas and elegantly confirmed the link between hop rates 
and the features of the potential landscape. Since this 
first quant itive test of Kramers' theory, a number of ex- 
periments have been published investigating how the hop 
rates of stochastic oscillations behave under modulation 
of the potentials jBHTO]. By asymmetrically modulating 
the depth of the potential wells and hence the phase be- 
tween the forcing and the activated transitions, Simon 
and Libchaber [6 observed synchronization of stochas- 
tic oscillations with a period equal to the mean Kramers 
time. 

On this micron scale, both fluctuations and friction 
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display long range hydrodynamic correlations that in- 
evitably couple the dynamics of suspended objects. For 
example, two colloidal particles in independent optical 
traps display correlations in their Brownian motions [11 . 
Whether the same hydrodynamic interaction could result 
in some degree of correlation in the thermally activated 
dynamics of nearby bistable systems is still an open and 
important question. Understanding how hydrodynamics 
can bias activated hops over energy barriers impacts the 
study of the slow dynamics of structural rearrangements 
in colloidal glasses and gels. 

In this Letter we use two colloidal spheres trapped in 
neighboring but separate bistable optical landscapes to 
reveal a significant difference between symmetric and an- 
tisymmetric simultaneous hops. Whilst we measure a sig- 
nificant biasing, the total number of simultaneous hops 
transpires only to be that expected from two random 
oscillators. This is due to the stochastic nature of our 
oscillators, whose dynamics are driven by thermal fluc- 
tuations alone. We show that the difference between the 
number of symmetric and antisymmetric hops reduces as 
the inverse separation between the two systems, show- 
ing that the strength of the phenomenon scales with first 
order in the hydrodynamic coupling. 

Sculpting of the bistable optical landscape is achieved 
by a holographic phase element generated by a gratings 
and lenses algorithm ^2] [13]. Holograms are imaged 
at the back aperture of a Carl Zeiss xlOO 1.3 NA oil- 
immersion microscope objective and focused into a mul- 
tispot array located at the focal plane (Fig. [T]). Each 
bistable potential is constructed by having two individ- 
ual focal spots placed in close proximity, such that a 
surmountable barrier forms between the two traps. We 
create two bistable potentials separated by a distance 5, 
which is also controlled holographically. Care is taken 
to maintain symmetric bistable potentials by adjusting 
the optical intensities of each trap. Stilgoe et al. have 
recently given a report on the complexity of creating op- 
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tical potentials using two optical traps [T4] . Our samples 
are prepared with Silica spheres (of radius, a = 0.4/im, 
Bangs Laboratories) in de-ionized, distilled water with 
a 1:10^ ratio. When preparing our samples we add a 
weak saline solution just before sealing so that any free 
debris is immobilised on the cover slip due to the re- 
duced Debye length [TF whilst still allowing enough time 
to capture the two spheres required for the experiment. 
Sample cells are constructed with a single concave mi- 
croscope slide (concave depth = 500 fim) and a square 
cover slip (thickness number 1.5), sealed with UV-curing 
optical adhesive. Using a few mW of frequency-doubled 
NdiYag laser light (532 nm), the trap stiffness in each 
trap is set at ~ IpN / iim~^ . A Prosilica GE650 camera, 
operating at a frame rate of 1.55 kHz is used to track 
the position of the spheres in the x — y plane. Using a 
center of mass algorithm, these in-plane displacements 
are measured with nm precision [16 . Time series of the 
spheres' center of mass, x{t) and y{t)^ are extracted in 
blocks of 10^ frames and accumulated for a total time 
period of 2-3 hours, significantly longer than the average 
between hops, which is of order one second. 
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FIG. 1: Holographic optical tweezers. A spatial light modu- 
lator is used to generate multiple traps at the focal plane of 
the xlOO oil-immersion microscope objective. Silica spheres 
are trapped with a frequency-doubled Nd:Yag laser, 5 /xm 
above the cover slip. Left inset shows the hologram used to 
landscape two side by side bistable potentials. A schematic 
of the two optically trapped spheres in independent bistable 
potentials is inset on the right. 

Since the sampling rate is much greater than the corner 
frequency in the power spectrum of particle fiuctuations, 
fc ~ 20 Hz^ the normalized histogram of particle posi- 
tions p(x, y) provides a direct measurement of the Boltz- 
mann distribution, p{x,y) = ex.-p[—U{x^y)/kBT]^ 
where ks-, T and Z are respectively the Boltzmann con- 
stant, the temperature of the surrounding fluid and the 
normalisation constant. By inverting p(x, y) we can di- 
rectly access the underlying optical potential in units of 
kBT (Fig.|and|3](b, d)). 

The obtained potentials are shown in Figure [5] (a - c) 
for three separations, s = 1.2, 2.0 and 2.8 /im. The po- 
tentials along the hopping axis for system a (blue o) and 
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FIG. 2: Optical potential landscapes experienced by two 800 
nm Silica spheres. Potential energy distributions are recov- 
ered from particle trajectories, x{t) and y{t), as described in 
the text. The data is binned into bins of size 10 nrri^ for a-b 
and 2 nm for d-f. 



system b (red o), presented in Fig. [2] (d-f), confirm that 
the optical landscapes of each system are not significantly 
perturbed as they approach. It is worth noting that, al- 
though the optical landscape is made of two optical traps 
whose centers are separated by 800 nm, the separation 
between the minima of the obtained bistable potential is 
only d ~ 400 nm. This is due to the sphere experiencing 
restoring forces from both optical traps, leading to an 
offset in the equilibrium position in each individual trap 
and shortening the hopping distance. 

Figure[3](a, c) shows a 10 5 sample trace of particle po- 
sition on the X-axis (the axis along which hopping occurs) 
taken at large separation, s = 3.6 /im. The correspond- 
ing potentials are given in Fig. |3](b, d). A single particle 
power spectrum is shown in Fig. [3] (e) for the shortest 
separation, s = 1.2 jim. Two separated timescales are 
visible, a low frequency component describing the hop- 
ping dynamics and a high frequency component describ- 
ing the fiuctuations around the local minima. For an 
isolated hopping system, the Fourier amplitude, xiv)^ of 
the particle's x position will be distributed with a double 
Lorentzian spectrum. 
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The first term is the standard spectrum for a bead of 
mobility mo = l/{67rr]a) (where r] = ImPas is taken 
as the viscosity) in an optical trap of strength /c, and 
thus, a corner frequency given by fc = mok/27T. The 
second term describes hopping events over a distance d 
and a jump probability per unit time given by fh- The 
hopping frequency can be anticipated using Kramers the- 
ory fh = cjo exp[— A/7//cb^] where AU is the height of 
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FIG. 3: (a, c) 10 s sample of the spheres' position along the 
X-axis (black line) for each system separated by s = 3.6 /urn. 
The block filtered (blue line) and digitised traces (red line) 
are also shown, (b, d) Corresponding potentials (black o) 
with fitted polynomials (red line), (e) Single particle power 
spectrum for s = 1.2 /jm along with a double Lorentzian fit 
as discussed in the text, (f) Correlation spectra for three 
separations evidencing hydrodynamic correlations in the high 
frequency region. 



the potential barrier and the prefactor ujq = rriy/KiK^ 
depends on the particle mobility mo and on the curva- 
tures of the underlying potential at the stationary points 
f^i.s = \d'^U /dx^\^^^ . Knowing the temperature, vis- 
cosity and particle radius means that all remaining pa- 
rameters can be directly evaluated from the features of 
the stationary points in the obtained potential energy. 
The resulting curve is shown as a solid line in Fig. [3] (e) 
and fits the measured power spectrum well. Searching 
for hydrodynamic correlations we plot in Fig. [s] (f) cor- 
relation spectra {x^{u)xi){h')) normalized by the average 
single particle power spectra (1/2) Xla^a 6(^a(^)^a(^))- 
A marked positive correlation at high frequencies gets 
larger as the two systems approach. This is the well 
known phenomenon by which the mobility of rigid mo- 
tions is higher than that for relative motions so that at 
short times the particles tend to move collectively pUj . 
However, at low frequencies dominated by hopping dy- 
namics, no sign of correlation is observable even when 
the two systems are close together. 

Despite the lack of any obvious sign of collective mo- 
tion in the low frequency correlations spectrum (Fig. |3] 
(f)), further analysis of our data reveals a significant cor- 
relation. Rather than averaging the correlation over all 
points, we examine only the instances where both par- 
ticles hop at the same time, either in the same direc- 
tion (symmetric) or opposite direction (antisymmetric). 
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TABLE I: Characteristic dwell times for each individual po- 
tential obtained by exponential fitting to the probability dis- 
tributions of the states dwell times, along with the mean value 
of raw dwell times and from Kramers formula. All numbers 
are in units of s. 



To examine these instances we first average the position 
time series in blocks of 60 frames (Fig. [3](a, c) blue line). 
This acts as a low pass filter with a cutoff frequency of 
about 25 Hz^ just higher than the trap corner frequency, 
filtering out the fast dynamics of fluctuations during a 
hop. The blocked data is then digitized into a two state 
variable x{t) G —1, +1, denoting the left and right states 
of the bistable potential (Fig. |3](a, c) red line). For each 
system we then extract a time series of hop events from 
— 1^+1 and +1 —1. This data represents a cumu- 
lative acquisition of 215 minutes at 1550 Hz^ with ~ 10^ 
transition events. The distributions of intervals between 
hops are fitted with an exponential decay, yielding the 
characteristic dwell time for each state. These measured 
dwell times are consistent with Kramers formula applied 
to the detailed shape of the underlying potential. Each 
bistable potential is fitted with an eighth order polyno- 
mial (Fig. [3] (b, d) red line), the second derivative of 
which gives the Kramers prefactor term ujq when evalu- 
ated at the stationary points. Numerical results of cj^^ 
are given in the third column of Table [Tj The dwell times 
are measured from x(t), the mean of which are given in 
the second column. The probability distribution of the 
dwell times are fitted with an exponential decay giving 
the first column in Table HI 

Each transition event is assigned a duration of one 
blocked sequence of frames, St = 39 ms. We count 
the number of occurrences of possible coincident events, 
namely p and n, where we define a p coincident as hav- 
ing occurred when both spheres hop in the same direction 
and an n coincidence when both spheres hop in opposite 
directions. Figure [4] (a) shows this nomenclature. It is 
expected that two uncoupled stochastic oscillators will 
show a finite number of coincidences, (A^), depending on 
the observed time period, T, and the temporal window, 
St^ in which we define a coincidence and the transition 
rates of the two oscillators, r~^, given by. 



{N) = 5tT{Tan) 



(2) 



For two uncoupled stochastic oscillators we expect an 
equal probability for symmetric and antisymmetric si- 
multaneous hops {p = n). Figure [4] (a) shows the mea- 
sured probabilities of p and n coincidences occurring over 
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FIG. 4: (a) Probability of coincidences over a range of system 
separations, s. p and n coincidences are defined by the inset 
cartoon. Coincidence events are counted as described in the 
text and presented as red o for p and blue o for n. Data is 
normalised by p + n. The solid lines represent the splitting 
in p and n assuming a linear dependence on hydrodynamic 
coupling with a strength of 3a /4s. (b) The total counted 
coincidences at each system separation normalised by the ex- 
pected number, (N), (Eq. [5]) with an error of iby^p + n/ (N). 



a range of system separations. With the two systems 
placed in close proximity to each other {s ^ 3a) the 
probability of observing a p coincidence is significantly 



greater than observing a n coincidence {p — n> {N) ~ 2 ) . 
At increasing separations, both p and n coincidence prob- 
abilities approach 0.5 ruling out the possibility that the 
synchronisation could be ascribed to asymmetries in the 
optical landscapes. We note that the total number of ob- 
served coincidences, (p+n)/(A^), is unchanged with sepa- 
ration (Fig. [4](b)). Coming back to Kramers' formula for 
single particle hopping rate, we notice that particle mo- 
bility appears as a prefactor to an energetic activation 
term. We may then conjecture that the rates of sym- 
metric hops and antisymmetric hops have a similar form 
with an activation term and a prefactor that would also 
depend respectively on the collective and relative mobil- 
ities pmiQ] . Moreover, the splitting between p and n 
seems to decay with distance as is the case for hydrody- 
namic splitting between collective and relative mobilities 
as shown by the 0.5 ± (3/4)(a/s) lines in Fig. [I] (a). 

In conclusion, the thermally activated jumps of two 800 
nm silica spheres in neighboring bistable optical land- 
scapes are shown to be coupled via hydrodynamic in- 
teractions. Due to the higher mobility of collective mo- 
tions, when two systems are in close proximity, there is 
a higher probability of observing a symmetric hop whilst 
antisymmetric hops are less common. We argue that the 
experimental environment studied here provides an ide- 
alised representation of interacting stochastic oscillations 
that occur in nature. It will be interesting to extend the 
present study to a larger ensemble of bistable systems. 
For example, looking for the emergence of more com- 
plex, cooperatively rearranging regions could aid the un- 
derstanding of the role of hydrodynamic interactions in 
the glassy dynamics of concentrated colloidal suspensions 
i20j. 
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